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ABSTRACT 

We report on a new approach to the calculation of thermodynamic functions for 
crossing-invariant models solvable by Bethe Ansatz. In the case of the XXZ Heisen- 
berg chain we derive, for arbitrary values of the anysotropy, a single non-linear 
integral equation from which the free energy can be exactly calculated. These equa- 
tions are shown to be equivalent to an infinite set of algebraic equations of Bethe 
type which provide alternatively the thermodinamics. The high-temperature ex- 
pansion follows in a sistematic and relatively simple way from our non-linear in- 
tegral equations. For low temperatures we obtain the correct central charge and 
predict the analytic structure of the full expansion around T = 0. Furthermore, we 
derive a single non-linear integral equation describing the finite-size ground-state 
energy of the Sine-Gordon quantum field theory. 
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The computation of thermodynamical functions for integrable models started 
with the seminal works of C.N. Yang and CP. Yang [1], of M. Takahashi [2] and 
M. Gaudin [3]. In Refs. [2,3] the free energy of the Heisenberg spin chain is written 
in terms of the solution of an infinite set of coupled nonlinear integral equations, 
derived on the basis of the so-called "string hypothesis" . In ref. [4] , we proposed a 
simpler way to solve the thermodynamics by means of a single, rigorously derived, 
nonlinear integral equation. This method can be applied to a wide class of models 
solvable by Bethe Ansatz (see ref. [5]). 

The XXZ hamiltonian for a periodic chain with 2L sites (we assume periodic 
boundary conditions but generalization to other b.c. is possible) takes the form 

2L 

HXXZ = -JY1 K^n+l + (^>l+l - COS7 («+i + 1)] (1) 

n=l 

The anisotropy parameter 7 is assumed here to lay in the real interval (0, vr) 
characteristic of the gapless regime. Hxxz is related to the diagonal-to-diagonal 
transfer matrix Ti{6) of the six- vertex model by 

Tl{0) '=° 1 - TTj^Hxxz + 0{9^) (2) 

2 J sm7 



where 



TUO) = R12R34 ■ ■ ■ R2L-1 2LR23R45 ...R2LI (3) 



and 



sin(7 — 9) , sin ^ 
a = , b = , c = 1 



(4) 
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6 is known as spectral parameter while a, b, and c are the conventional six-vertex 
Boltzmann weights. 



The free energy per site is defined as usual 



(3 L^oo 2L ■ 
and from eq. (2) we read 



/(/3) = -- lim —log [Tr fe-f'H^-A] (5) 






N 

Tl{2(3/N) 



where (3 = (3Jsm'~f. Hence we can write 

fiP) = -i hm ^ Jim \ogZLN{2P/N) (6) 

where Zin(9) = Tt[T]:^{9)] is the six-vertex partition function on a periodic 
diagonal lattice with L x N sites. The two limits in eq. (6) cannot be interchanged 
since the degeneracy oiTi{0) = 1, that is 2^^, is strongly L— dependent. However, 
the crossing invariance of the six-vertex i?— matrix implies that under a rotation by 
7r/2 of the entire lattice plus the substitution 6 ^ '-/ — 6, the numerical value of the 
partition function does not change (see e.g. ref.6). Therefore Zin{9) = Z]s[l{i — 0) 
and 

/(/3) = -]- hm ^ lim \ogZNdl - 2(3/N) (7) 

where Zj^^iiO) = Tr \T^{6)] . Now, as the well-known Bethe Ansatz solution tells 
us, T/v(7) has a non-degenerate largest eigenvalue for any finite A^. Then the two 
limits in eq. (7) commute and one finds [7] 



-2(3f{(3)= lim logAr"(7-2/?/Ar) 



where K^°'^{d) denotes the largest eigenvalue of Tn{6). Using the results of ref. [8] 
for the eigenvalues of the diagonal-to- diagonal transfer matrix, the free energy 



becomes 



where 






EN{P)-2Nlog 



sm{2p/N) 
sin 7 



(9) 



N 



e(A, P) = 0(A, + 27/2, 7/2 - P/N) - 0(A,- - 27/2, 7/2 - ,3/iV) 



,., , ., sinhfix + A) 

0(A, x) = 2 log — — — 

^\mi\ix — X) 



(</.(0,x) = 0) 



(10) 



The real numbers Ai, . . . , Aat are the roots of the Bethe Ansatz equations 



iV 



0(A,-,/3/iV) + 0(A,-,7-/3/iV)] =Y.^{^^ - Afc,7) +vr(2j -iV- 1) (11) 



fe=i 



associated to the largest eigenvalue Ajy°'^(7 — 2j3/N) (we assumed A^ even). We 
define as usual the counting function [9] 



zn{\) 



27r 



1 ^ 
0(A, ^/A^) + 0(A, 7 - /9/A^) - ^ E <^(^ - ^fc' ^) 



A:=l 



'121 



which by construction enjoys the properties -2iv(A) = —zjq{—\) = Z]^{X) and 



zn{\j)=J-{N + 1)/2 



(13) 



Our aim is to derive an integral equation for zn{X). Let us start by assuming that 
the root distribution becomes the continuus density o"c(A) (normalized to 1) in the 
limit A^ -^ 00. Then it could be identified with the derivative of zpf{\) and at 
the same time used to replace the discrete sum in eq. (12) with an integration. 



For large but finite N we whould therefore obtain the following approximate linear 
integral equation for o"c(A) 

+00 
27rac(A) =0'(A,/3/iV) + 0'(A,7-/9/iV)- I d^i<p'{\- ^x,^)a,{^i) (14) 



Fourier transforming leads to the explicit solution 

sinh(7rA/7) 



o"c(A) = -2c(A) , Zc{\) = — arctan 



_sm{7ip/-fN)_ 



(15) 



We see that o"c(A) is strongly peaked at A = for large N, reflecting the singularity 
which develops at the origin in the N ^ 00 limit of the source term in eqs. (11), 
(14). The density picture correctly describes the BA roots wherever o"c(A) is of 
order 1. From eq. (15) we then find as validity interval 



|A|<0(V^) (16) 

which shrinks to zero when N ^ 00. Roots |Aj| > 0{y(3/N) have a spacing of 
order larger than 0{1/N) and cannot be descried by densities. In particular, the 
roots with largest magnitudes have finite N —^ 00 limits spread by 0(1) intervals 
(we checked this fact numerically too). Therefore, contrary to the usual situation 
[9], we must go beyond the density description to obtain a bulk quantity like the 
free energy per site. The sum in the last term of eq. (12) can be written as a contour 
integral 

' '^' ^(A-^,i.)^^^5!l±^M£)i (17) 



{^}i:^(^-^..7)=/ 



27riiV " dz 

3^' 



where ajv = exp(27riA^ZAr) and the contour in eq.(17) encircles counterclockwise all 
the real roots A^, A; = 1, ..., A^. Notice that ajv(Afc) + 1 = 0, A; = 1, ..., iV. Inserting 



eq.(17) in eq.(12) and integrating by parts yields 

+00 
\ogaN{\))=tN (f){\,P/N) + (l){\,^-P/N) +t / — 0'(A-x,7)Imlog[l+aN(x+iO)]| 



TT 

— CO 

(18) 
Convoluting eq.(18) with the kernel 

K-' = 6{X) - p{X) 

where 

+00 

— 00 

we find 

+00 

aiv(A) f l + aiv(/i + iO) 

log^-— = / rf/ip A-/i log- — . 20 

—00 
which is the sought nonlinear integral equation. (Here a^ = exp{2iTiNzc)). Notice 



that ajv(A) = l/ajv(A). We evaluate the sum in eq. (10) by a similar procedure, 
with the result 

+00 

T (f^\--[? (a\ T?(f^\ ^ /"^\ sinh7rA/7 cos7r/^/7iV l + a]v(A + iO) 

Ljv(/3) = En[P) - Ec[l3) = - dX — — -— ^ -, ,, log ■ 



7 J cosh 27rA/7 - cos 2v:j3/'^N 1 + aN{X + iO) 

—00 

(21) 
where E(.{(3) is that part of the energy that follows from the root density (15), that 
is 

+00 +00 
^ , ^, /",,/. X /.x ,. /" (i/c sinhfvr — 7)A; sinh(7 — 2/5/A^)A; , , 
Ec{(3) = / rfAe(A,7)a,(A) = -2N / ^ . '' , V , ^ ^ 22) 

J J k smh7rA;cosh7A; 

—00 —00 

If eqs. (20) and (21) are analytically continued to the axis ImA = ±7/2 they 
assume the same structure of those derived in ref. [12] by different methods. 



Inserting eqs. (21) and (22) in eq. (9) yields 



fm = Exxz{l)+p-'HP) 



(23) 



where Exxzil) is the ground-state energy of (1), namely 



Exxzil) = 2J 



+CX) 



COS 7 — sin 'J dk 



sinh(7r — '~f)k 
sinh7rA;cosh7A; 



(24) 



while L[f3) is the N ^ oo limit of (21), that is 



+C30 



L{(3) = Im d\ 



log[l + a(A + iO)] 
7sinh7r(A + iO)/7 



(25) 



Similarly, the function a(A) in this expression is the A^ ^ oo limit of ajv(A). It is 
therefore the solution of the N ^ oo limit of eq. (20), that is 



-iloga(A) 



2ttI3 



+O0 



7sinh7rA/7 



2 / rf/ip(A-/i)Imlog[l + a(/i + iO)] (26) 



Thus the calculation of the free energy (5) has been reduced to the problem of 
solving the single nonlinear integral equation (26) and then evaluating the integral 

(25). 

We find from the N = oo limit of eq.(18) an alternative form for our eq. (26) 



+00 



27r^ 



z{\) = /3q{\) + —^Im / (iyU0'(A — /i, 7) logcos7rz(/x + iO) (27) 



where 



g(A) = - 



1 / sinh2A 



TV \cosh 2A — cos 27 



— coth A 



(2^ 



and z{\) = 2:^loga(A) is related to the original counting function 

z{\) = lim N[zNi\)-hign\] (A real) (29) 

By the residue theorem we then obtain 

z{X) = hW - ^ E [<^(^ - ^^•' ^) - '^(^' ^)] (30) 

jeTL 

where the real numbers ^j, defined by z{^j) = j — 1/2, j G TL, can be identified with 
the N -^ oo limit of the original BA roots Ai, A2, • • • , A^v-The new algebraic BA 
equations z{C,j) = j — 1/2 embody all the information about the XXZ thermody- 
namics. Eq. (30) shows that z{X) has periodicity iir and has, as unique singularity 
on the real axis, a simple pole at the origin with residue — /S/tt. The new algebraic 
BA equations z{^j) = j — 1/2, j G TL embody all the information about the XXZ 
thermodynamics. They form an infinite set of Bethe Ansatz type algebraic equa- 
tions equivalent to our non-linear integral equations (26) and (27) . We want to 
stress that these roots are discrete although we have already set iV = 00. The 
root spacing is actually of order T. In the T — ^ limit we recover the continu- 
ous distribution of roots associated with the antiferromagnetic ground-state [9]. In 
terms of the roots ^fc, the free energy takes the form 

f{(3) = Jcos7 - /3"' ln2 + ^ 5^ log ([1 - cos 27]/ [cosh 2^, - cos 27]) (31) 

The roots ^j have an accumulation point at .^ = 0. That is, ^j = — /3/[7r(j — 1/2] 
for j -^ 00. For low temperatures, the largest root, .^1, is of order (7/vr) lnf3. 

Let us now study /(/3) for high temperatures. When (3 is small it is convenient 
to use the alternative form (27) plus the uniform expansion 



z{\) = Y,P'hW (32) 

fc=i 



Since the residue of z{X) is linear in j3 and z{X) is an odd function, we have 



61(A) ^~°-^ + 0(A), 6fc(0) = ik>2) (33) 



Inserting eq. (32) in eq. (27) we find, up to fourth order in (3 



61(A) = g(A) 62(A) = ^0"(A,7) 



47r 

a- „ 

sm 7 



.3(A) = HX) = ^ i - -i- ^"(A.7) + ^r'(A.7) 



(34) 



Then, from eqs. (23) and (25), 

f{f3) = _/3-iln2 + Jcos7-/5j2(l + icos2 7)+/5^^^cos7 + 0(/53) (35) 

which indeed agrees with the high T expansion [2]. We want to remark that eq. 
(27) generates the functions bk{X) recursively, with easy integrations which involve 
only delta functions and derivatives thereof. It is indeed a very efficent way to 
recover the high-temperature expansion from the Bethe Ansatz solution [5] . 

We shall consider now the low temperature regime. When /5 ^ 1 eq. (26) 
indicates that z{X) ~ (3, so that log[l + a(A)] ~ 27riz(A), at least as long as 
|A| < (7/vr) log/3. At dominant [3^1 order eq. (26) linearizes in the same way as it 
does for small p. Inserting then z{X) ^ Pq{X) in eq. (25), yields hm^_+tx) L{P) = 0. 
Therefore eq. (23) tells us that 



f{(3)^='^Exxz{l) + 0{r') (36) 



The contributions of order P ^ and smaller come from values of A larger than 
(7/77) log /9, where the previous assumption log a ~ 0{f3) does not hold anymore. 
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It is then convenient to introduce the new function 



7 AttP 
A{x) = a(A) , X = \ log 

TT 7 



Then eqs. (25) and (26) reduce, in the /5 — > oo hmit, to 



+ 00 

L{p) = X f dx e-^^/^Im log[l + A{ 

■Kj3 J 



X + ii 



(37) 



(3^ 



and 



+CX) 



-i\ogA{x) = -e-^""/^ + 2 I dyp{x-y)Im\og[l+A{y + iO)] (39) 

— oo 

Fortunately, it is not necessary to solve this integral equation to calculate the 
integral (38). In fact, we can use the following lemma. 

Lemma. Assume that F{x) satisfies the nonlinear integral equation 



X2 

—i log F{x) = (p{x) + / dyp{x — y) Im log[l + F{y + i 



0)] (40) 



Xl 



where ^{x) is real for real x and xi, X2 are real numbers. Then the following 
relation holds 



Im / dx^\x) log[l + F{x + iO)] = ilm [ip{x2) log(l + Fa) - v?(a;i) log(l + Fi)] 

Xl 

+ iRe [£(Fi) - £(F2)] 
where Fi^2 = F(a;i,2) and £ is a dilogarithm function 



(41) 



i{t) = du 





log(l + u) \ogu 
u 1 + u 



(42) 
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To prove the lemma we consider the relation 
i{F2)-i{Fi)= f dx Uog[l + F{x)]-^logF{x) - logF{x)^\og[l + F{x)] 



and then use eq. (40) and its derivative to substitute log-F(a;) and dlogF{x)/dx. 
After taking real part and a little algebra this yields eq. (41) (related identities 
were used in ref. [12]). 

The integral in eq. (38) may now be exactly calculated by invoking the lemma 
with F{x) = A{x), ip = — exp(— 7ra;/7) and xi = — cxd, X2 = +cxd. We have 
A{xi) = 0, A{x2) = 1 and 



+ CX) 



IT O 



Then the free energy for low temperature reads 

/(/?) = Exxzil) - TTT—r^ + o{r^) (43) 

bJ sm7 

in perfect agreement with refs. [2,10,11]. As we have just shown, both the high and 
the low temperature leading behaviors of the free energy can be derived without 
much effort from our non-linear integral equation (26). The higher order correc- 
tions for high T can be obtained in a very sistematic way [5] . The situation for the 
o(/9~^) terms in the low T expansion (43) is more involved. Qualitatively, however, 
it is rather easy to establish, from the asymptotic behaviour of the kernel p{X — ^) 
in eq. (26) and from the in periodicity implied by eq. (30), that these higher order 
terms must be integer powers of T'^"' and T^iy'^^^) . 

Let us now consider the problem of calculating the (properly subtracted) 
ground state energy E{L) of 2D integrable massive models of Quantum Field 
Theory in a finite 1-volume L. This subject recently received much attention 

12 



in connection with Perturbed Conformal Field Theory. In this contest E{L) is 
known as ground-state scahng function. Our starting point is the lattice regular- 
ization of such models provided by the light-cone approach to integrable vertex 
models [8]. Specializing to the six-vertex model, the ground state energy, on a ring 
of length L formed by A^ sites, takes a form similar to eq. (21): 



+CX) 

ENiL)-EciL) = ^lm f d\ 

— 00 



sech ^(A - 0) - sech ^(A + 6) \og[l+aN{\+ii 



where 



Ar2 
Ec{L) = — 



+ 00 



-27r + / d\ 



0(A + 20,7/2) 



7 cosh 7rA/7 
and ajv(A + iO) obeys an equation like (20) with the new source term 



(44) 



(45) 



logac(A) = 2arctan 



sinh 7rA/7 
cosh TT 6/7 



(46) 



The L— dependence in these expressions is hidden in the light-cone parameter B, 
which tends to infinity, in the continuum limit iV — ^ 00, as [8] 



e = - log ^ 



(47) 



with m, the physical mass scale, held fixed. Define now 



e(A) = — hm loga7v(7A/7r) 

iV-^00 



(4J 



then from eqs. (20), (42-45) we find 



+CX) 



E{L) = lim [^iv(^) - Ec{L)] = -m I —sinh Aim log 



l + e 



-e{\+iO) 



(49) 
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and 

+00 

e(A) = mL sinh A + 2 f dfiGiX-fx) Im log [l + e~<''+'^^] (50) 

—00 

where G{X) = {'~f/7i)p{'~fX/7r). Thanks to the correspondence between the hght- 
cone six-vertex model and the sine-Gordon (or Massive Thirring) model [8], eqs. 
(49) and (50) give the ground state scaling function E{L) of the sine-Gordon field 
theory on a ring of length L.(Here m stands for the fermion mass). This is a 
rigorous and simpler alternative to the standard Thermodynamic Bethe Ansatz 
[14]. Notice that all UV divergent terms are contained in Ec{L) (eft. eq. (45)). 
E(.{L) also contains the finite, scaling part m^Lcoi^ 7r^/27 (also found in ref. [13] 
by completely different means), which gives the ground-state energy density in the 
L -^ 00 limit. 

Let us now study eqs. (49) and (50) for mL ^ 1 and mL <^ 1. For small 
values of m,L the calculation closely parallels that for low T in the thermodynamics 
of the XXZ chain. Shifting A ^ A — log m,L and applying the lemma leads to 

E{L)"'^='-^ + 0{l) (51) 

showing the expected value c = 1 for the UV central charge. The large mL regime 
easily follows from eq. (50) by iteration. One finds the behavior typical of a massive 
quantum field theory 

E(L) "^^^^ -'^KiimL) + 0(e-2-^) (52) 

71 

with Ki{z) the modified Bessel function of order 1. 

We would like to remark that, contrary to the traditional thermal Bethe Ansatz 
[2,3], our approach does not relay on the string hypothesis on the structure of 
the solutions of the BA equations characteristic of the Heisenberg chain. This 
makes our approach definetly simpler. Most notably, the whole construction of 
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the thermodynamics no longer depends on whether ^j-n is a rational or not, unlike 
Takahashi approach [2]. This applys equally well to the problem of the ground 
state scaling function of the sine-Gordon field theory, since the standard TBA 
approach requires the string hypothesis. 

Generalizations of the present TBA approach to higher spin chains as well as to 
magnetic chains and quantum field theories associated to Lie algebras other than 
A\ (nested BA solutions) are relatively straightforward provided crossing symmetry 
holds [5]. 
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